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Abstract 


The theory of solar gamma-ray "Ine production is reviewed and new 
calculations of line production yields are presented* Observations, carried 
out with gamma-ray spectrometers on OSO-7, HEAO-1, HEAO-3 and SMM are reviewed 
and compared with theory. These observations provide direct evidence for 
nuclear reactions in flares and r'urnish unique information on particle 
acceleration and flare mechanisms. 
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I. INTRODUCTION 

Nuclear reactions In solar f lures take place between flare accelerated 
protons and nuclei and the ambient solar atmosphere. Several reviews of the 
early work on this subject are available (Dolan and Fazio 1965, Lingenfelter 
and Ramaty 1967, Cheng 1972). Using reasonably accurate and complete nuclear 
data, Lingenfelter and Ramaty (1967) have carried out the first detailed 
calculation of the expected nuclear reaction rates in flares and predicted 
observable fluxes at Earth of the products of these reactions, gamma-ray 
lines, neutrons and nuclear fragments in the solar energetic particles, 

Nuclear gamma rays from solar flares were first observed by Chupp et al. 
(1973) with a Nal spectrometer flown on board the seventh Solar Orbiting 
Observatory (0$0-7),« Gamma-ray lines at 0.51 MeV, 2.22 MeV, 4.44 MeV and 6.13 
MeV were observed from the August 4, 1972 flare. The lines at 0.51 and 2.22 
MeV were also seen during the decay phase of the August 7, 1972 flare. These 
two flares were among the largest ever observed. 

A considerable amount of theoretical and interpretative work has been 
carried out on the August, 1972 bservations (Ramaty and Lingenfelter 1973, 
Reppin et a’l. 1973, Wang and Ramaty 1974, Forrest, Chupp and Suri 1975, 

Ramaty, Kozlovsky and Lingenfelter 1975, ^^ang 1975, Kanbach et al. 1975, Chupp 
1976, Lin and Hudson 1976, Bai and Ramaty 1976, Ramaty and Crannell 1976, 
Crannell et al. 1976, Kozlovsky and Ramaty 1977, Ramaty, Kozlovsky and Suri 
1977 Ibragimov and Kocharov 1977, Lin and Ramaty 1978, Crannell, Crannell and 
Ramaty 1979, Ramaty 1979, Ramaty et al. 1980). These studies, together with 
additional gamma-ray data (Chupp, Forrest and Suri 1975, Chupp 1976, Suri et 
al. 1975), and energetic-particle (Kohl, Bostrom and Williams 1973, Webber et 
al. 1975), hard X-ray (van Seek, Hoyng and Stevens 1973) and microwave (Croom 
and Harris 1973) data have lead to the following broad outline of the origin 
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and implications of solar ganina rays. 

Gamma-ray lines from fiOlar flares result from the interaction of flare 
accelerated protons and nuclei with the solar atmosphere. The emission 
measure of flare heated material, i.e. the density squared of the hot gas 
times its volume which is determined from X-ray observations, is insufficient 
to produce any measurable amount of nuclear burning. The accelerated 
particles, on the other hand, produce neutrons, positrons, k mesons, 
radioactive nuclei and excited nuclear levels, whose captures, annihilations, 
decays, and deexcitations lead to observable gsmma-ray lines. 

The strongest line in solar flares is at 2.223 MeV from neutron capture 
on hydrogen. The neutrons are produced mainly from the disintegration of ^He 
and heavier nuclei and occasionally in proton-proton collisions resulting in 
high-energy neutrons and tt mesons. High-energy neutrons could be directly 
detected near Earth (Lingenfelter et al= 1965, Lingenfelter and Ramaty 1967) 
and preliminary reports of their observation have just become available (Chupp 
and Forrest 1981). 

The site of the nuclear reactions in the solar atmosphere is as yet 
unknown. Nevertheless, calculations (Wang and Ramaty 1974) indicate that the 
bulk of neutrons with initial velocity vectors pointing towards the 
photosphere are thermal i zed in the photosphere and subsequently captured on 
either or ^He. Capture on produces the 2.223 MeV gamma-ray line, but 
capture on ^He results in tritium without emitting photons. The removal time 
of thermal neutrons from the photosphere, on the order of a minute, can be 
measured by comparing the time dependence of the intensity of the 2.223 MeV 
line to that of a prompt nuclear line (see below). This removal time depends 
on the photospheric ^He abundance, as does the capture probability on which 
determines the flux of the 2.223 MeV line. Gamma-ray line observations. 





therefore, can measure the abundance of %e in the photosphere. Because the 
2.223 Me‘- line is formed at a larger depth in the solar atmosphere than the 
prompt gamma-ray lines, for flares close to the limb of the Sun, the 2,223 MeV 
line is substantially attenuated relative to the prompt emissions. 

A variety of prompt gamma-ray lines are produced from nuclear 

deexcitations. The most important discrete lines (Ramaty, Kozlovsky and 

Lingenfelter H»75, 1979) are at 6.129 MeV from 4.438 MeV from 2.313 

MeV from 1.779 MeV from , 1,634 MeV from 20Ne, 1.369 MeV from 

1.238 MeV and 0.847 MeV from ®^Fe, all produced primarily by direct excitation 

of these nuclei, and at two lines, 0.478 MeV from \i and at 0.431 MeV from 
^Be, which result from fusion reactions, He(a,p) Li and He(a,n) Be . 

role of these fusion reactions for producing gamma-ray lines in astrophysics 
was first pointed out by Kozlovsky and Ramaty (1974). Nuclear deexcitations 
also produce Doppler broadened lines which together with many unresolved lines 
produce a significant gamma-ray continuum, in particular In the 4 to 7 MeV 
region (Ramaty, Ko;ilovsky and Suri 1977, Ibri5gimov and Kocharov 1977), 

Because of the short lifetimes of most excited nuclear levels, nuclear 
deexcitation radiation is an excellent tracer of the nuclear interaction rate 
of energetic particles in solar flares. This rate is directly proportional to 
the instantaneous number of accelerated particles in the interaction region, 
which, in turn, is determined by the acceleration mechanism and the losses 
suffered by the particles. Through line shapes and Doppler shifts, prompt 
nuclear deexcitation lines also give information on the geometry of the 
interacting energetic particle beam (Ramaty and Crannell 1976, Kozlovsky and 
Ramaty 1977). 

The interaction models of energetic protons and nuclei in solar flares 
can be crudely classified as thin- and thick -target models (e.g. Ramaty, 
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Kozlovsky and Llngenfelter 1975). In the thin-target model the nuclear 
reactions are produced by energetic particles which escape from the 
Interaction region at the Sun. These particles can be detected In the 
Interplanetary medium. Furthermore, If sufficient thin-target nuclear 
reactions take place, their fragmentation products should also be 
detectable. On the other hand, In the thick-target model the nuclear 
reactions are produced by particles as they slow down In the solar 
atmosphere. These particles, and their spallation products, become 
thermal Ized and mixed back Into the solar atmosphere. Gamma-ray lines and 
neutrons can, nevertheless, be seen from thick-target Interactions. 

The ratio of the flux In a prompt nuclear component (e.g. the 4.44 MeV 
line) to the flux in the 2.223 MeV line depends on the Interaction model and 
on the energy spectrum of the accelerated particles. For a given energetic 
particle spectrum, this ratio Is larger for the thin- target than for the 
thick-target model. Using this result and Information on the energy spectrum 
derived from interplanetary particle observations (Van Hollebeke, MaSung and 
McDonald 1975), Ramaty (1979) suggested that the gamma rays from the August 4, 
1972 flare were produced predominantly in thick -target interactions. From the 
absolute fluxes of the gamma ray lines, it is possible to deduce the energy 
deposited by the accelerated nuclei in the solar atmosphere. For the August 

4, 1972 flare, the energy deposited by the protons and nuclei amounts to 

/ 

several percent of the energy depo^-ted by the electrons which make the 
Impulsive hard X-rays (Lin and Hudson 1976). Nevertheless, protons could 
deposit their energy In regions which are not accessible to electrons because 
they have a longer stopping range in the ambient medium than the < 100 keV 
electrons. 

Positrons In solar flares result from the decay of ti'*’ mesons and various 
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radioactive nuclei produced by the nuclear reactions* The half lives of the 
iinportant positron emitters range from values less than 1 second to over 20 
minutes and they produce positrons of energies from about. 0,1 MeV to several 
tens of MeV. After their production, the positrons are decelerated to 
energies less than several hundred eV where they annihilate. The deceleration 
Is due to Interactions with the ambient solar atmosphere, and hence the 
deceleration time depends on the density and magnetic field of the medium In 
which the positrons annihilate. The positrons can annihilate with free 
electrons to produce two 0.511 MeV gamma rays per annihilation, or they may 
form a positronlum atom. This atom Is similar to the hydrogen atom except 
that the proton Is replaced by a positron. Positronlum atoms also annihilate 
Into gamma rays: 25% of the annihilations are from the singlet spin state 

producing two 0,511 MeV photons, and 75% of them are from the triplet state 
which annillates Into three photons of energies less than 0.511 MeV,^, Triplet 
positronlum can annihilate before It Is broken up by collisions only if the 
density of the ambient medium 1s less than about 10^® cm“^. Observation of 
the characteristic 3-photon positronlum continuum, therefore, would provide 
information on the density of the annihilation site (Crannell et al, 1976). 
Information on the temperature of this site would be obtained from the 
measurement of the width of the 0.511 MeV line. If the positrons annihilate 
below the transition layer, I.e. at temperatures less than lO^K. the width the 
0.511 MeV line should be less than 3.5 keV. 

Following the OSO-7 observations, solar gamma-ray lines were seen with 
the Nal spectrometer on the first High Energy Astrophysical Observatory (HEAO- 
1) (Hudson et al. 1980), with the Nal spectrometer on the Solar Maximum 
Mission (SMM) (Chupp et al. 1981, Chupp 1981), and with the 6e spectrometer on 
HEAO-3 (Prince et al. 1982). The two HEAO detectors, designed to detect 
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cosmic gamma rays, observed solar gamma-ray lines through their shields and 
thus were sensitive to only the strongest one or two lines. The HEAD 
detectors are no longer operational. The SMM spectrometer has already seen 
lino amission from many solar flares and thus demonstrated that nuclear 
reactions of flare accelerated particles take place commonly in solar 
flares. At the time of this writing (early 1982) the detector is still 
operational and should continue to observe for several more years. The ten 
solar flares, from which gamma-ray lines were definitely seen so far, are 
discussed in the present paper. 

In addition to gamma-ray lines and neutrons, nuclear reactions also 
produce energetic nuclear fragments (e.g, 3)^^ 3^e, Li, Be, B). The first 

attempt to measure the %e abundance in energetic solar particles was made by 
Schaeffer and Zahringer (1962). By using mass spectroscopy of material from 
the Discoverer 17 satellite, thes? authors found a ^He/^He ratio of -- 0.2 for 
the November 12,1960 flare. Subsequent measurements (Hsfeh and Simpson 1970, 
Anglin, Dietrich and Simpson 1973a, Dietrich 1973, Garrard, Stone and Vogt 
1973) have revealed the existence of a class of solar particle events in which 
the ^He/^He ratio is substantially larger than in the ambient solar 
atmosphere. 

Enrichments of ^He in energetic particle populations (for example, the 
galactic cosmic rays) have been generally attributed to nuclear reactions 
between the energetic particles and the ambient medium. But, as first pointed 
out by Garrard, Stone and Vogt (1973), this interpretation of the solar ^He 
enrik'hments, in its simplest form, is inconsistent with much of the ^He 
data. If the ^He enrichments are due to nuclear reactions of the energetic 
particles, then they should be accompanied by similar enrichments in and, 
to a lesser degree, in ^H. Such enrichments, however, are not observed. 
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Several schemes have been proposed to overcome this difficulty. These 
rely on the kinematics and angular distributions of the reaction products 
which favor the preferential escape of ^He (Ramaty and Kozlovsky 1974, 

Rothw^n 1976) and the thermonuclear destruction of and In a model In 
which the energetic products of the nuclear reactions are confined to thin 
filaments and Interact with each other (Colgate, Audouze and Fowler 1977), 

But, as proposed by Fisk (1978), the enhanced ^He abundance In solar energetic 
particles could be due to preferential heating and acceleration of ambient 
atmospheric ^He. The observation of energetic ^He, therefore, cannot be used 
as indication for nuclear reactions unless the ^He is accompanied by at least 
some other fragmentation product. No convincing observations of such products 
have yet been reported. 

Nuclear reactions of accelerated particles could cause modifications of 
solar surface isotopic abundances, Thus, from the analysis of lunar surface 
material, Kerridge (1975) found a secular increase of the solar wind 
ratio, and Fireman, De Felice and D'Amico (1976) reported a measurable 
abundance on the lunar surface which they believe should be due to 
implantation by the solar wind. Even though solar surface nuclear reactions 
could, in principle, produce these isotopes, it Is unlikely that this has 
indeed happened, because the necessary nuclear reaction rates on the ancient 
Sun would have to exceed the present rate (determined by the gamma-ray 
observations) by many orders of magnitude (Kerridge et al. 1977). 

Thus, the only convincing evidence to-date for nuclear reactions in the 
solar atmosphere are the gamma-ray line and neutron observations. The 
observed Tine energies and line ratios are fully consistent with nuclear 
reactions produced by particles of energies in excess of several MeV. As 
already mentioned, thermonuclear burning makes no measurable contribution to 
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these reactions. In Section II we review the Interaction models for the 
production of gamma rays In energetic particle reactions and present new 
calculations of the production rates. In Section III we review the 
observational data on gamma-ray lines, we compare them with theory and we 
discuss their Implications, Because the neutron observations are still very 
preliminary, we defer their analysis to future studies. We summarize our 
conclusion In Section IV. 

II, NUCLEAR REACTIONS IN SOLAR FLARES 

In this section we first summarize the formalism of the thin- and thick- 
target interaction models used in ^ic’lculating nuclear reaction yields in solar 
flares, and we discuss the energetic particle spectra and compositions that we 
use in these calculations. We then evaluate, for the various models, the 
neutron and 2.223 MeV photon productions, the positron and 0.511 MeV photon 
productions and the productions of the various prompt nuclear deexcitation 
lines. 

A. Interaction Models and Properties of the Energetic Particles 

We consider first the thin-target model. Here nuclear reactions take 
place between accelerated particles and a cold ambient medium In an 
interaction volume from which particles escape with negligible energy loss and 
with an escape probability which is energy independent and the same for all 
types of particles. Let Nj(E,t)dE be the instantaneous number of energetic 
particles of type j In the volume having energies per nucleon in dE around E, 
and the density of ambient particles of type 1. The Instantaneous reaction 
rate between accelerated particles of type j and ambient particles of type i 
is given by 


u 


q^j(t) » /* Nj(E,t)$(E)o^j(E)dE . (1) 

where a(E) is particle velocity and is the energy dependent cross 

section of the reaction considered. 

In the thin-target model the energetic particles that escape from the 
interaction region can* in principle, be observed in the interplanetary medium 
by detectors on spacecraft. The following relationship exists between the 
interplanetary particles, Nes(^j(E), and the instantaneous number of particles 
in the thin-target interaction volume; 

H^s^^j(E) » T‘^ / dt Nj(E,t) dt . (2) 

Here T is the escape time from the volume and the time Integral is over the 
duration of the nuclear interactions. By integrating equation (1) over time, 
by substituting equation (2), and by summing over all 1 and j that contribute 
to a particular reaction product (e.g. neutrons, positrons, excited levels), 
we obtain the time- integrated nuclear reaction yield of that product in the 
thin- target model; 

As opposed to the thin-target model, in which the particles escape from 
the Interaction region and can be observed in space, the thick-target model is 
one in which the particles produce nuclear reactions as they slow down In the 
solar atmosphere. In this model, the ambient density in the interaction 
region is expected to be quite high, (i.e. the region could be close to or 
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even Inside the photosphere). Since particles, in general, are not 
accelerated in high density regions* it is reasonable to assume that in the 
thick-target model the acceleration takes place outside the interaction 
volume. Thus, let fJj(E)dE be the time-integrated number of particles of type 
J with E in dE incident on the interaction region. The time-integrated 
nuclear reaction yield can then be written as 



(4) 


where mp is the mass of the proton and (dE/dx)j is the energy loss rate per 
unit path length (measured in (MeV/nucleon)/(g/cm2) ) of accelerated particle j 
in the ambient solar atmosphere. In the present calculations we take (dE/dx)j 
equal to the energy loss rates of charged particles in a neutral medium, 



[ 1 + 


"He "He , 

"h > i 


(5) 


Here {dE/dx)j^H energy loss rates of particle j in H 

and He, respectively, and m^g is the mass of ^He. For the abundances of Table 
1, the term in the square brackets is approximately 1.13 and essentially 
independent of energy. Using the tabulations of Barkas and Berger (1964) and 
Northcliffe and Schilling (1970), (dE/dX)j^^^ can be approximated by 

(^) “ (Zfiff/^^j 630 E”°‘®(MeV/nucleon)/(g/cm^ 

j»H 


( 6 ) 
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where (Pierce and Blann 1968) 

Zgff = Z[ 1- exp (-137p/Z^''^] , (7) 

and Z and A are the atomic and mass number of particle j. For computational 
purposes, It Is convenient to Invert the order of Integration In equation (4), 


0 = 


"’p'^^j^^i/Hh) ^ (dE/dx)J^ 


/ dE' R. 
E ^ 


(E') 


( 8 ) 


A detailed discussion of solar energetic particle spectra and 
compositions based on Interplanetary observation are given In ..her article 
in this volume (Forman, Ramaty and Zwelbel 1982). Various forms of 
accelerated particle spectra have been used In previous treatments of nuclear 
reactions In solar flares (Lingenfelter and Ramaty 1967, Ramaty, Kozlovsky and 
Lingenfelter 1975, R*?.maty 1979). These are power laws In kinetic energy. 


N 




“S 


esc.J<E) 0^ Nj (E» = t 3 .-s 

j C 


E>E. 


E<E 

- C 


exponentials In rigidity. 


(9) 


Nesc,j(E) or Nj(E) = exp( - R./R^) dRj/dE 


( 10 ) 
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and Bessel functions, 


KgjC,j(E> <»• [ 2<3p/(mpC“'''>>'^^! • <“> 

In these expressions the Bj's are proportional to the abundances of energetic 

2 _ 

particles j, Rj = (A/Z)jp 1s particle rigidity, p = /E(E + 2mpC ) Is particle 
momentum per nucleon and K 2 is the modified Bessel function of order 2 (e.g. 
Abramovltz and Stegun 1966). The parameters s and for power laws, Rq for 
the exponential In rigidity, and aT for the Bessel function, characterize the 
spectrum of the energetic particles. 

Recent observations of energetic protons and a-particles on spacecraft 
(McGuire, von Rosenvinge and McDonald 1981) Indicate that power laws cannot 
fit the observed energy spectra of these particles. This result Is consistent 
with earlier studies (Freier and Webber 1963) which have shown that 
exponentials In rigidity provide a better fit to the solar proton data than do 
power laws In kinetic energy. Therefore, we shall present new calculations of 
nuclear reaction yields only for the exponential spectrum of equation (10) and 
the Bessel function spectrum of equation (11). Calculations for power laws 
can be found In previous papers (Ramaty, Kozlovsky and Lingenfelter 1975, 
Ramaty, Kozlovsky and Surl 1977). 

In Figure (1) we show the energy spectra of protons and a-particles 
observed by detectors on the IMP 8 spacecraft from the June 7,1980 flare (R. 
McGuire, private communication 1981). This flare Is one of the bes*", studied 





gamma-ray flares observed by the SMM spectrometer (Chupp et al. 1981). As can 

1 /2 

be seen, the curves In Figure 1, giv^n by dj/dE « e K2C2(3p/(m cuT)) ' ] 

H 

with aT s 0.015, provide a very good fit to the data. The June 7, 1980 flare 
is the only gamma-ray flare observed to-date for which good interplanetary 
energetic particle spectra have been reported. But even for this flare it is 
not entirely clear that all the interplanetary particles shown in Figure 1 
were produced by the same flare as the one which produced the gamma-ray lines 
because of evidence for multiple injection of particles from the Sun (T. Von 
Rosenvinge, private communication 1981). 

Enuation (li) was shown (M. Lee private communication 1978, Ramaty 1979, 
Forman, Ramaty and Zweibel 1982) to be the solution of a Fokker-Pl anck 
equation for stochastic Fermi acceleration with acceleration efficiency 
coefficient, a, and escape time, T, which are independent of particle energy 
and particle charge. While there is no guarantee that such a simple 
acceleration model is appropriate for solar flares, we feel that the spectrum 
of equation (11) can be adequately used for the analysis of the presently 
available gamma-ray data, particularly since this spectrum also provides a 
good fit to the observed interplanetary particle energy spectra, as can be 
seen from Figure 1 and from the more detailed analysis of McGuire, von 
Rosenvinge and McDonald (1981). 

The particle abundances that we use in the present calculation are given 
in Table 1. The ambient medium abundances (Cameron 1981) are listed in the 
second column. From interplanetary observations, it is well known that the 
energetic particle abundances vary substantially from one flare to another 
(e.g. Forman, Ramaty and Zweibel 1982). In the present calculations we use 
two sets of energetic particle abundances. The first set, denoted by EPl, is 
given in the third column of Table 1, and the second set, EP2, is identical to 
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the ambient medium abundances of the second column* The energetic particles 
EPl are significantly richer In heavy elements than EP2. 

B. Neutron and 2.223 MeV Photon Productio n 

The strongest line observed In nearly all gamma ray flares Is that at 
2.223 MeV from neutron capture on hydrogen, ^H(n,Y)^H» Several theoretical 
studies have been made of neutron production In solar flares (Lingenfel ter et 
al. 1965, Lingenfel ter and Ramaty 1967, Ramaty, Kozlovsky and Lingenfel ter 
1975) and of 2.223 MeV photon production from the capture of these neutrons In 
the solar atmosphere (Wang and Ramaty 1974, Kanbach et al. 1975). The neutron 
production cross sections have been discussed In considerable detail by 
Ramaty, Kozlovsky and Lingenfelter (1975). These cross sections have recently 
been updated (B. Kozlovsky and R. E. Lingenfelter, private communication 1981) 
with the addition of many new reactions that Involve all the Isotopes listed 
in Table 1. These new cross sections, to be published elsewhere, are used in 
the present calculations. 

The calculated neutron production yields, Q^, In the thin and thick- 
target models are shown In Figure 2. These calculations are normalized such 
that the number of escaping protons of energies greater than 30 MeV, 

Noe/. r»(>30 MeV), and the number of protons Incident on the thick target above 

f P 

the same energy, (>30 MeV), are both equal to unity. The energetic 
particles abundances are given by EPl (see Table 1). For relatively flat 
energetic particle spectra, corresponding to large values of Rq or aT, the 
bulk of the neutrons are produced In reactions between protons and 
a-particles. For steep particle spectra, given by the Bessel function at 
small aT, the large neutron yields result from reactions between a-particles 
and heavy nuclei. These large neutron yields are absent for the rigidity 
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spectra because particles with Z>2 have lower energies per nucleon at the 
same rigidity and hence produce less nuclear reactions. These effects were 
discussed in more detail by Ramaty, Kozlovsky and Lingenfelter (1975) in 
connection with the camparison of neutron production by particles with spectra 
that were either power laws in energy or exponentials in rigidity. 

In Figure 3 we show partial neutron production rates in the thick target 
model for energetic particle spectra given by equation (11) and abundances 
again given by EPl. The curves labelled p, a, CNO and Ne-Fe represent, 
respectively, neutron yields of energetic protons, a-part1cles, C, N and 0 
nuclei, and nuclei from Ne through Fe interacting with the ambient solar 
atmosphere. As can be seen, except for the very steep particle spectra (small 
otT) the neutrons result mostly from a-part1cles and protons. For very flat 
spectra, the neutron yield of protons includes an important contribution from 
the reaction p+p > p+n +tt . 

In Figure 4 we show the energy in accelerated particles required to 
produce one neutron in the thick- target model. Here W is defined by 


W = lA. / dE' E' R.(E') , (12) 

j ^ E J 

where fl.(E) Is given by equation (11) and the abundances are given by EPl. As 
J 

can be seen from Figure 4, for aT = 0.015, appropriate for the June 7, 1980 
flare (see Figure 1), more than Vz ol' the accelerated particle energy Is 
contained In particles of energies greater than 1 MeV/nucleon and this 
fraction Increases with increasing «T. A much lower fraction of the energy 
content, however, resides in particles of higher energies. 
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The values of W/Q„ given in Figure 4 should be considered as lower limits 
because the particle energy loss rate (equation 6) Is valid only for slowing 
down of test particles in a neutral medium. In a fully ionized medium, the 
si owing-down rate of test particles is larger by about a factor of 3(e,g. 
Ramaty 1979). Furthermore, collective effects, such as the generation of an 
induced magnetic field by a beam of particles (Hoyng, Brown and Van Beek 1976, 
Colgate 1978) would increase the energy required for the production of a given 
amount of neutrons. 

The 2*223 MeV gamma-ray line is formed by neutron capture on in the 
photosphere. To study this line formation, Wang and Ramaty (1974) have 
carried out a detailed Monte-Carlo simulation in which a distribution of 
neutrons was released above the photosphere, and the path of each neutron 
after its release was followed. For isotropic neutron release, any initially 
upward moving neutron escapes from the Sun. Some of the downward moving 
neutrons can also escape after being backscattered elastically by ambient 
protons, but most of these neutrons either are captured or decay at the Sun. 
Because the probability for elastic scattering is much larger than the capture 
probability, the majority of the neutrons are thermal ized before they get 
captured. Since the thermal speed in the photosphere (where most of the 
captures take place) is very much smaller than the speed of light, the energy 
of the gamma rays is almost exactly 2.223 MeV (Taylor, Neff and King 1967), 
and the Doppler-broadened width of the line is very small (<100 eV) (Ramaty, 
Kozlovsky and Lingenfelter 1975). The energy of the line observed with the 
high resolution Ge detector on board HEAO-3 from the November 9, 1979 flare is 
2.225 ± 0.002 MeV (Prince et al. 1982). 

The bulk of neutrons at the Sun are captured either on or on ^He. 
Whereas capture on yields a 2.223 MeV photon, capture on ^He proceeds via 
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the radiationless reaction ^He(n,p)% and hence produces no photons. The 
cross sections for capture on and ^He are 2,2xl0~^® cm^ and 3.7xl0~26 
respectively, where ca^, is the velocity of the neutron (for details 
see Wang and Ramaty 1974). Thus, If the ^He/H ratio In the photosphere Is 
5x1 comparable to that observed In the solar wind (Geiss and Reeves 1972) 
and In the chromosphere (Hall 1975), nearly equal numbers of neutrons are 
captured on ^He as on H. 

The results of the Monte-Carlo calculations of Wang and Ramaty (1974) are 
presented in Figure 5. In these calculations an Isotropic distribution of 
monoenergetic neutrons of energy is released above the photosphere. The 
solid lines are the probabilities for the various Indicated processes. As can 
be seen, the capture and loss probabilities Increase with Increasing energy, 
because higher energy neutrons penetrate deeper into the photosphere. This 
reduces their escape probability and leads to a shorter capture time, thereby 
reducing the decay probability. The probability for loss on ^He almost equals 
the capture probability on protons. The escape probability is greater than 
0.5, because all Initially upward moving neutrons were assumed to escape from 
the Sun. Note that the sum of all probabilities equals 1. 

The dashed curves In Figure 5, are neutron-to-photon conversion 
coefficients evaluated (Wang and Ramaty 1974) for specific emission angles e 
between the Earth-Sun line and the vertical to the photosphere and given 
neutron energies. At low neutron energies and 9 near zero, f2,2 1^ close to 
the capture probability on protons. This means that gamma rays from low- 
energy neutrons observed close to the vertical escape essentially unattenuated 
from the Sun. At higher energies and at larger angles, however, there is 
significant attenuation of the gamma rays due to Compton scattering in the 
photosphere. Therefore, for flares close to the limb of the Sun, the 2.223 
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MeV line should be strongly attenuated in comparison with other nuclear 
deexcitation lines which are likely to be produced at higher altitudes In the 
solar atmosphere than the 2.223 MeV line. As we shall see in Section III A, 
this limb darkening Is clearly seen in the SMM data. 

The time Integrated flux, or fluence, of 2.223 MeV photons at Earth 
resulting from neutron capture at the Sun can be written as 

4»(2.223 MeV) - 

where d = 1 AU and neutron-to-2.223 MeV photon conversion 

coefficient averaged over the neutron energy spectrum. For 6=0, Ramaty, 
Kozlovsky and Lingenfelter (1975) find that 2 from about 0.1 to 

0.14, depending on the neutron energy spectrum. However, In the thick-target 
model, the neutron angular distribution Is probably not Isotropic and the 
neutrons could be produced In the photosphere (Kanbach et al. 1975). These 
effects should Increase by as much as a factor of 2. In our subsequent 
discussion we denote by the conversion coefficient for vertical escape by 
Fg 2 ^®^ • calculations now definitely justify new and more accurate 

calculations of neutron capture in the photosphere. These have not yet been 
carried out. Nevertheless, using the presently available calculations, vje 
believe that «(0) should not exceed 0.3. 

C. Positron and 0.511 MeV Photon Production 

The 0.511 MeV gamma-ray line resulting from positron annihilation has 
been observed from several solar flares. Positrons In flares are produced in 
energetic particles Interaction with the ambient solar atmosphere. A number 
of theoretical studies have been made of positron production In such 
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Interactions (Llngenfelter and Ramaty 1967 » Ramaty, Kozlovsky and Lingenfelter 
1975) and the positron slowing down and annihilation (Crannell et a1, 1976, 
Bussard, Ramaty and Drachman 1979). In the present paper we give the results 
of new calculations (B. Kozlovsky and R. E. Llngenfelter, private 
communication 1981) of positron production based on a large number of b’*' 
emitters produced In nuclear reactions that involve all the Isotopes listed In 
Table 1. The results are given In Figure 6, where we show the ratio Q^/Q^ for 
the thin- and thick-target models. Here Q+ is calculated from equations (3) 
and (8), respectively, with energetic particles given by equation (11). The 
positron yields shown In this figure represent total yields. Because of the 
finite half-lives of the various emitters, however, In a short observation 
time of a transient event, fewer positrons than indicated In Figure 6 are 
available for 0.511 MeV line production. This effect is shown in Figure 7, 
where dQ+/dt Is the Instantaneous production rate of positrons from a burst 
(6-funct1on in time) of 3'*’-emitters produced at t = 0. 

The positrons produced In nuclear reactions have Initial energies ranging 
from about 0.1 MeV to tens of MeV, depending on the production mode. In a 
thick-target, these positrons are rapidly slowed down to energies of tens of 
eV where the majority of them form positronium atoms. The slowing down time, 

tg, plus the positronium formation time, tpos, are shown in Figure 8 as a 
function of Initial positron energy for an ambient medium of temperature lO^K, 
degree of Ionization n and ambient density lO^l cm"^ and 10l3 cm“3 (from 
calculations of J. M. McKinley, private communication 1981), The dependence 
of (tg + tpQg)"^ on density is essentially linear. 

The annihilation of positronium atoms is quite rapid (Heltler 1954). If 
formed In the singlet state { 25 % of the time), positronium annihilates at a 
rate of 8x10^ sec"^ into two 0.511 MeV gamma rays. In the triplet state 
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(formed 75% of the time) It annihilates at a rate of 7x10® sec“l j^to three 
photons of energies less than 0.511 MeV. This 3-photon continuum Is, In 

principle, observable If the density of the ambient medium Is less than *-10^® 
cm"^. At a higher density collisions break up triplet positronlum before It 
can annihilate (Crannell et al. 1976), 

We define a pos1tron-to-0.511 MeV photon conversion coefficient, ?q 
analogous to g discussed above, such that the time Integrated flux, or 
fluence, of 0.511 MeV photons at Earth resulting from positron emitter 
production at the Sun is given by 

<j.(0.511 MeV) = ^ 

Various effects Influence the value of gj^ . If all the emitters 
produced In the flare decay during the observation period and If all the 
resultant positrons annihilate In this period, then ?q g^ ranges from 0.5 to 
2, depending on the fraction of positrons that annihilate via positronlum. 

But If the observation period is shorter than the decay halflives of the 
dominant p'*' emitters, or If some of the positrons escape from the Sun into low 
density regions where their annihilation time is long, 7 q gj^ can be 
substantially lower than the above values. For limb flares, there may be 
significant Compton scattering of the 0.511 MeV photons in the photosphere if 
the nuclear reactions. themselves take place in the photosphere. This will 
also lower the value of Tq 5^ • 

Another observable of considerable Interest is the width of the 0.511 MeV 
line. The dependence of this width on temperature, degree of Ionization and 
density was studied In considerable detail by Crannell et al. (1976) and 
Bussard, Ramaty and Drachman (1979). Using their results, we find that the 
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O, 511 MeV line should be narrower than about 3.5 keV If the temperature of the 
annihilation site Is less than 10®K. At higher temperatures the width should 
vary as with a full width at half maximum of about 11 keV at 10®K. 

P , Prompt Deexcitation line Production 

A variety of gamma-ray lines are produced In solar flares from de- 
excitation of nuclear levels. Figure 9 shows the spectrum from these de- 
excitations, calculated (Ramaty, Kozlovsky and Lingenfelter 1979) by employing 
a Monte-Carlo simulation for an energetic particle population Interacting with 
an ambient medium. The energetic particle spectrum Is proportional to 
E“2 (equation 9 with Eg « 0) and both the ambient medium and the energetic 
particles at the same E have a photospheric composition. The shapes of the 
lines are evaluated by taking Into account nuclear kinematics and data on the 
differential cross-sections of the reactions. The results of the simulations 
are binned Into energy Intervals ranging from 2 to 5 keV (as Indicated in the 
figure), consistent with the resolution of a Ge gamma-ray spectrometer. 

Two line components can be distinguished in Figure 9. A narrow component 
resulting from the deexcitation of ambient, heavy nuclei excited by energetic 
protons and a-part1cles, and a broad component from the deexcitation of 
energetic heavy nuclei Interacting with ambient H and He, The nuclei 
responsible for the emission of the strongest narrow lines are Indicated In 
the figure. 

In addition to these strong narrow lines, there are many other weaker 
narrow lines, which together vf*th the broad component produced by heavy 
accelerated particles, merge Into the underlying continuum. Above '•* 4 MeV 
most of the radiation Is from C, N, and 0, while below about 3 MeV the 


principal contributors are Mg, SI and Fe, It has been shown (Ramaty, 

Kozlovsky and Surl 1977, Ibragimov and Kocharov 1977) that the bulk of the 
gamma-ray flux observed between 4 and 7 MeV from the August 4, 1972 flare was 
of such nuclear origin rather than electron bremsstrahlung or other continuum 
emission processes. In addition, a substantial fraction of the photons In the 
1 to 2 MeV band could he nuclear radiation resulting from an enhanced 
abundance of Ne, Mg, SI and Fe in the energetic particles (Ramaty et al. 1980 
and Section III A). The 4 to 7 MeV energy band, referred to In the SMM 
observations as the "main channel window", can provide a direct and sensitive 
measure of the Interaction rate of protons and nuclei In solar flares (Section 
III A, B). 

Because the cross sections for excitation of nuclear levels have 
different energy dependences from those of neutron and positron production, 
the calculated ratios of nuclear deexcitation line yields to the neutron yield 
depend strongly on the assumed spectra of the accelerated particles and on the 
Interaction model. As an example, in Figure 10 we show the ratio Q(4.44)/Qp 
for the two Interaction models as a function of aT for the Bessel function 
spectrum of equation (11). Here Q(4.44) is the yield of the narrow 4.44 MeV 
line calculated from equations (3) or (8), for thin or thick targets, 
respectively, with energetic particle abundances given by EPl (see Table 1). 

As can be seen, Q(4.44)/Qf, generally decreases with increasing energetic 
particle spectral hardness, reflecting an increased neutron production and 
decreased 4.44 MeV photon production by particles of high energies. Likewise, 
this ratio Is lower for thick targets than for thin targets (except for very 
steep particle spectra) because the energy losses harden the particle spectra 
In the thick target. For very steep (small «T) Bessel function spectra, the 
thick-target ratio exceeds that for thin targets. This results from the 
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effect of the energy losses in the thick-target which suppress the heavy 
particle fluxes relative to the proton and a-particle fluxes. At low particle 
energies, the heavy particles contribute significantly to neutron production 
but not to the production of narrow 4*44 MeV photons. 

As already mentioned, the photon energy band from 4 to 7 MeV is an 
important measure of the Interaction rate of protons and nuclei in solar 
flares. In Figure 11 we show the ratio of the photon yield in this band, Q(4- 
7), to the neutron yield Q^. Here Q{4-7) is calculated from equations (3) ur 
(8) for thin or thick- targets, respectively, together with a Monte-Carlo 
simulation similar to that employed in the evaluation of the gamma-ray 
spectrum of Figure 9. In this calculation we use both the EPl and EP2 
energetic particle abundances. 

The same trends as In Figure 10 are also evident for Q(4-7)/Q^ in Figure 
11. Q(4-7)/Qp decreases with increasing aT and is larger in the thick-target 
than in the thin-target model. The variation of Q(4-7)/Q^ with energetic 
particle composition, however, is not very large, because several particle 
species contribute simultaneously to both Q{4-7) and Q^. 

The ratio Q(4.44)/Q(4-7) is shown in Figure 12 for thin and thick targets 
and the energetic particle abundances EPl and EP2. As can be seen, this ratio 
does not depend strongly on energetic particle spectrum because of the similar 
energy dependences of all prompt line production cross sections. However, 
Q(4.44)/Q(4-7) is smaller for the thin target than for the thick target, 
because for the latter the contribution of the heavy particles is suppressed 
by the energy losses. Likewise Q{4.44)/Q(4-7) is lower for EPl than for EP2 
because the former contains more heavy particles relative to prsjtons than does 
the latter. 

The shapes and peak energies of prompt nuclear gamma-ray lines are 
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sensitive to anisotropies in the energetic particle angular distributions. 
Thus, Ramaty and Crannell (1976) have evaluated the shift in the peak of the 
6.129 MeV lines resulting from energetic particle beaming, while Ramaty, 
Kozlovsky and Lingenfelter (1979) discusseu the splitting of 4.438 MeV line 
that is observed (Kolata, Auble and Galonsky 1967) when the excited nuclei 
are produced by a proton beam perpendicular to the direction of observation. 

A unique test of energetic particle beaming was proposed by Kozlovsky and 
RfimSity (1977). This concerns the ^Be* and \i* deexcitation lines a* 0.431 
and 0.478 MeV produced in the reactions ^He(«,n)^Be* and '^He(a,p)^Li*, 
respectively. Here the stars indicate nuclei in excited states. The shapes 
of these lines are shown in Figure 13. As can be seen, if the angular 
distribution of the energetic a-particles is isotropic, the Doppler broadening 
of the two lines is so large that they blend into a single feature that 
cannot, in general, be observed in the presence of a strong continuum. If, 
however, the a-particles are beamed, the line widths are much less than in the 
Isotropic case and two discrete lines can be seen. In particular, if the 
direction of observations is perpendicular to the beam (as in Figure 13) the 
lines appear close to the rest energies f 0.431 MeV and 0.478 MeV. 

The fluence at Earth in the 4-7 MeV channel is given by 

<t>(4-7 MeV) = Q(4-7 MeV) ?4.7/(47rd^) , (15) 

where y is the conversion coefficient from nuclear deexcitations to 
photons observed. If the excited nuclei are produced isotropically and if 
there is no attenuation of the photons, = 1 . However, if the protons 
and nuclei form a beam pointing away from the observer, then from the 
Monte-Carlo calculations described above we find that y « 0.8 . Thus, 
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provided that, there Is not much attenuation, we expect that 0.8 ij ^ i 1 , 
depending on the geometry of the Interacting particles. 

Ill, IMPLICATIONS OF GAMMA RAY OBSERVATIONS 

In this Section we consider the Implications of the gamma-ray line 
observations on the nature of the Interaction model {thin- or thick-target), 
on the accelerated particle spectrum, on the energy deposited by the particles 
and the number of particles that Interact to produce the gamma rays, on the 
timing of the acceleration, on the photospheric ^He abundance, and on the 
beaming of the energetic particles. 

The solar flares with observed gamma-ray lines are listed In Table 2. 

The August 4 and August 7, 1972 events were observed by Chupp et al. (1973) 
with the spectrometer on board OSO-7. The data given In Table 2 for the 
August 4, 1972 flare Is from Chupp (1976), except for (|»(4-7 MeV) which Is from 
the analysis of Ramaty, Kozlovsky and Surl (1977). We have multiplied the 
time averaged fluxes given In these references by 553 sec, the observation 
time of gamma rays from this flare (Chupp 1976). Because of Earth occultatlon 
of the orbiting gamma-ray detector, however, the total duration of gamma-ray 
emission from the August 4 flare was longer than this time Interval. From the 
analysis of Wang and Ramaty (1975), we estimate that all fluences In Table 2 
for the August 4 flare should be Increased by about a factor of 2, but without 
a significant modification of <|)(4-7)/(t>(2.22). 

Only the 2.22 and 0.51 MeV lines were observed from the August 7, 1972 
event (Chupp 1976) because the detector was behind the Earth during the flash 
phase of the flare. These observations clearly demonstrate the delayed nature 
of these two lines: at a time when all prompt emissions were very small, the 
2.22 and 0.51 MeV lines were still observable. 
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The data for the flares of July 11, 1978 and November 9, 1979 are, 
respectively, from Hudson et al. (1980) and (Prince et al. 1982). 

The data for the other flares listed In Table 2, still believed to be 
preliminary, were observed with the gamma-ray spectrometer on SMM (Chupp 
1982). Because of the very high intensity of the June 21, 1980 flare, the SMM 
detector saturated during the flash phase of this event. Delayed 2.22 and 
0.51 MeV lines, as well as high energy neutrons (Chupp and Forrest 1981), were 
observed from the June 21 flare. 

The ratios of the photon fluence in the 4 to 7 MeV channel to that in the 
2.22 MeV line, shown in Table 2, follow directly from the data except for the 
July 11, 1978 flare where <j>(4-7) is determined using a theoretical ratio 
Q(4.44)/Q(4-7) = 0.25. The fact that this ratio is model dependent (Figure 
12) leads to some uncertainty in the determination of (|)(4-7) for this flare. 

A. Interaction Model, Energetic Particle Spectrum, Number and Energy Content 

We first consider the flare of June 7, 1980 for which there are both 
gamma-ray line observations and interplanetary particle measurements. The 
combined analysis of these data imply that for this flare the bulk of the 
gamma-ray line emission results from thick target interactions. Furthermore, 
most of the energetic protons and nuclei that produce the gamma-ray lines 
remain trapped in the solar atmosphere and only a small fraction of them 
escapes into the interplanetary medium. 

The location of the June 7 flare at N12W74 indicates that it was well 
connected magnetically (e.g. Van Hollebeke, Ma Sung and McDonald 1975), so 
that particles escaping from the Sun could be observed in interplanetary space 
near the orbit of the Earth. Indeed, several observations of energetic 
particles have been reported (von Rosenvinge, Ramaty and Reames 1981, Evenson, 
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Meyer and Yanagita 1981, Pesses et al. 1981). Based on those, the number of 
protons of energies greater than 10 MeV released into the Interplanetary 
medium, Ngsc,p^^^® has been estimated (von Rosenvinge et a1. 1981) to be 
- 1031. As can be seen from Figure 1, the spectrum of these protons is well 
fit with the Bessel function of equation (11) with aT equal to 0.015 (R. 
McGuire, private communication 1981). This spectral form also fits the a- 
particle spectrum with essentially the same oiT. For aT = 0.015, Nesc,p^^^^ 
MeV) « 10^^ implies Nesc,p^^^^ ~ 5x10^^. If the gamma rays were produced 
by thin- target interactions, then from 'jquation 3, with the numerical results 
of Figure 2 and 2 ” (appropriate for a thin target), the observed 
2.22 MeV line fluence (Table 2) implies that n^T « 7.4x10^^ cm^sec"^. This is 
equivalent to a matter traversal for 30 MeV/nucleon particles of « 13g/cm^. 

The large abundances of spallation products (^H, ^H, Li, Be, B) that would 
result from such a long path length are not observed from solar flares (e.g. 
McGuire, von Rosenvinge and McDonald 1977). This indicates that the gamma-ray 
lines observed from the June 7,1980 flare were probably not produced in thin- 
target interactions. In the thick-target models, on the other hand, the 
spallation products that accompany the production of gamma-ray lines, are 
slowed down in the solar atmosphere and hence are not expected to be seen in 
the interplanetary medium. 

Analysis of the ratio of the 4 to 7 MeV fluence to the 2.223 MeV line 
fluence for the June 7, 1980 flare also suggests that the observed gamma rays 
were produced in thick-target interactions, not thin-target. This can be seen 

4 

as follows: 

From equations 13 and 15, with ^ = 1 , the neutron-to-2.223 MeV photon 
conversion coefficient, 2 ’ written as 





The numerator can be obtained from theory (see Figure 11) If aT Is known. The 
Interplanetary particle data suggests that oiT « 0.015 (Figure 1) provided that 
all the observed particles were Indeed produced in the gamma -ray flare and 
that the flare particle spectrum Is not greatly modified by the escape process 
and Interplanetary propagation. The denominator in equation (16) is from the 
gamma-ray data, <|>{4-7)/<}>(2.22) « 1.74 (Table 2). Then for thin target 
Interactions, Q(4-7)/Q„ a 0.7, hence ^2 2 ~ * while for thick target 

Interactions, Q(4-7)/Qn « 0.25 hence ^2 2 ~ calculations of 

Wang and Ramaty (1974) and Kanbach et al. (1981) we estimate that for the 
location of the June 7, 1980 flare, fg 2^^2 2'®^ " flare, 

therefore, gio) would have to be ~ 0.67, if the gamma rays were made In 
thin-target interactions, and 0.23 If they were made In thick-target ones. 

The former value Is clearly Inconsistent with the range of values 

of 2 ^®^ obtained from the Monte-Carlo simulations of neutron production In 

the solar atmosphere (Section II B). The latter value, for the thick target, 

Is quite consistent with these calculations provided the neutrons are produced 
In the photosphere and/or their initial angular distribution Is skewed 
downwards toward the photosphere. 

We proceed now to analyze the rest of the data listed In Table 2. We 
note the relatively small variability of (|)(4-7)/<ji(2.22) from one flare to 
another. The exception Is the limb flare of April 27, 1981 for which the 2.22 
MeV line Is strongly attenuated l[^ Compton scattering as the photons emerge 
from the photosphere at large angles to the local normal (Wang and Ramaty 
1974). We assume that In all of these flares, as In the June 7 flare, the r 

gamma rays are produced by thick-target interactions. This Is justified | 

i 





because different Interaction models for different flares would not be 
consistent with the relative constancy of <ji{4-7)/(|>(2.22) shown in Table 2. 
Furthermore, we assume the same neutron- to-2. 223 photon conversion coefficient 
for all flares, except for the correction due to flare location. We 
use ?2 2 ( 0 ) = 0 , 23 , a value consistent with both theory and the June 7 
particle and gamma-ray data. 

Table 3 lists the flares with available (j?{4-7)/t(2.22) ratios and the 
April 27, 1981 flare. The second column gives an estimate of fg 2/^2 2^®^ 
obtained from the flare locations and the calculations of Wang and Ramaty 
(1974) and Kanbach et al. (1981). The third column lists values of Q(4-7)/Q^ 
deduced from equation (16), while the fourth column provides the value of aT 
obtained from these ratios and Figure 11. We note the relatively small 
variability of aT from one flare to another, a direct consequence of the 
constancy of (j)(4-7)/<j>(2.22) . This result is consistent with the observations 
of McGuire, von Rosenvinge and McDonald (1981) who find a range of aT‘s for 
protons and a-particles in interplanetary space which essentially overlaps 
that deduced for the gamma-ray flares. Particle acceleration evidently 
produces energy spectra that do not vary much from flare to flare. A similar 
conclusion has been obtained by Van Hollebeke, Ma Sung and McDonald (1975). 

We cannot deduce the aT for the limb flare of April 27, 1981 because of 
the strong attenuation of the 2.223 MeV line by Compton scattering in the 
photosphere. The small variability of aT, however, allows us to assume an aT 
for this flare. We take aT = 0.019, equal to that for the August 4, 1972 
flare, because both ev.ents have similar durations of gamma-ray emission and 
approximately equal fluences in the 4 to 7 MeV channel. 

Using the aT's listed in column 4, we calculate, in the thick-target 
model, the number of particles that interact in the solar atmosphere and the 
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energy deposited by them. The results, based on the numerical values of 
Figures 2 and 4 are given In columns 5 and 6 of Table 3. Here, W(>1 
MeV/nucleon) Is the energy deposited by particles of energies greater than 1 
MeV/nucleon and Rp{>10 MeV) Is the number of protons of energies greater than 
10 MeV Incident on the thick-target. For the June 7, 1980 flare, the number 
of protons of energies greater than 10 MeV that interacts at the Sun exceeds 
that observed In Interplanetary space ('^lO^l) by about two orders of 
magnitude. In contrast, for the August 4, 1972 flare, the number of protons 
above 10 MeV observed in the interplanetary medium 3xlo35, Lin and Hudson 
1976) exceeds the number that interacts at the Sun (Table 3) by more than an 
order of magnitude. 

The very large number of interplanetary particles observed from the 
August 4, 1972 flare could have produced the observed gamma rays by 
thin-target Interactions, as proposed by Lin and Hudson (1976). These authors 
have assumed a flatter energetic particle spectrum than that given by aT = 
0.019 in Table 3. However, the interplanetary particle spectrum from the 
August 4 flare Is only very poorly known because several interplanetary shocks 
were present at that time. Rather than using the interplanetary observations, 
we would now argue that the relative constancy of (j)(4-7)/(j)(2.22) supports the 
same interaction model for all gamma-ray flares, and hence a thick-target 
model for the August 4, 1972 flare. This result is consistent with the 
average Interplanetary proton spectrum observed (Webber et al. 1975) from 
August 2 to 11, which can be well fitted with equation (11) with «T = 0.02, in 
good agreement with the value of oT given in Table 3 for a thick-target. 

The energy depositions of the protons and nuclei given in Table 3 range 
from about 5x10^8 erg to 2.5x10^® erg. However, as pointed out in Section II 
B, these should be considered as lower limits only, because of additional 
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energy losses to the Ionized component of the solar atmosphere and possible 
collective effects (e,g. Colgate 1978). Nevertheless, we can compare the 
energy deposition of the nucleonic component with that of >25 keV electrons 
deduced from hard X-ray observations In a nonthermal model (Lin and Hudson 
1976). The energy deposition of the electrons ranges from about 2 x 1Q29 erg 
for small flares to 10^^ ergs for the August 4, 1972 flare (Lin and Hudson 
1976). We see that the nucleonic component could be responsible for the 
deposition of at least several percent of the total flare energy. 

A final argument that supports the thick-target Interaction model comes 
from the analysis of the 0.511 MeV line from positron annihilation. In Table 
4 we list four flares for which there Is either an obi,^ervat1on of, or an upper 
limit on the fluence In this line. For the August 4, 1972 flare the data Is 
from Chupp (1976), for the dune 7 and July 1, i960 flares It Is from Chupp 
(1982) and for the April 27, 1981 flare it Is from D. Forrest (private 
communication 1981). Using the otT's of Table 3 and the results of Figures 6 
and 11, we calculate the ratios Q+/Q(4-7 MeV). These are shown in column 4 of 
Table 4. In column 5 of this table we give the values of (|)(0.51) calculated 
from the values of <{)(4-7 MeV) given in Table 2 and with ?q , the p'*' 
emitter-to-0.511 MeV photon conversion coefficient, a free parameter. The 
observed fluences or upper limits are given in column 6. By comparing the 
calculated and observed values of (|)(0.51), we see that for the August 4 event 
there is good agreement If 5], is about 0.7 which is consistent with the 
theoretical expectation discussed In Section II C. A similar value of gj^ 
could also account for the July 1, 1980 data. For the June 7, 1980 
flare, fn ci would have to be less than about 0.4, a value consistent with 
the short observation period (50 sec) which does not allow the complete decay 
of all the positron emitters produced by the nuclear reactions (see Figure 


i 



i».r* 


34 


7). To account for the observed upper limit, the value of for the April 

27, 1981 flare, however, must be less than about 0,.2 even though the 
observation period here (~ 2000 sec) Is sufficiently long for the decay of 
essentially all of the positron emitters. It appears that the best 
explanation of the absence of the 0.511 MeV line In the April 27,1981 limb 
flare Is Compton scattering In the photosphere (Ramaty, Lingenfelter and 
Kozlovsky 1982). Because the neutrons have relatively long stopping ranges, 

Compton scattering of the 2.22 MeV line In limb flares is expected whether or 
not the nuclear reactions take place In the photosphere. But because the 
positrons have generally shorter ranges, the 0.511 MeV line will be Compton 
scattered In such flares only if the nuclear reactions themselves take place 
In the photosphere. This result. If substantiated by further studies, should 
be a strong argument for the validity of the thick-target model. 

B. Time Dependences 

The time dependences of the gamma-ray lines contain Important Information 
on a variety of questions In solar physics. 

The time dependence of the strongest discrete line from flares, the 2.223 
MeV line from neutron capture, Is determined by the time history of the 
neutron production as well as by the removal rate of neutrons from the 
photosphere where they spend most of their time (--I minute) between production 
and radiative capture (Section II B). 

il 

The delay between neutron production and 2.223 MeV photon release has 

been unmistakably observed In several flares (Chupp et al. 1973, Hudson et al . I 
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1980, Chupp et al. 1981, Prince et al. 1982). Such an observation entails the ji 
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comparison of the time history of a prompt photon flux, for example that In 
the 4 to 7 MeV channel, with that of the 2.223 MeV line. This has been 
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carried out in detail for the Oune 7, 1980 flare (Chupp et al, 1981), 

According to Chupp (1982), the characteristic neutron removal time from the 
photosphere, as deduced from the gamma-ray data, Is on the order 50 seconds, 
consistent with theory (Section II B) and a photospheric ^He/H ratio of 
~5xl0"^. As we shall see below (Section IIIC), the fact that the delay 
between neutron production and capture Is not much shorter than 50 sec can 
place an upper limit on the photospheric ^He abundance. 

The time dependence of the 0.511 MeV line Is determined by the production 
rate of the a"*" emitters, the decay rate of these emitters (Figure 7) and the 
slowing down and annihilation time of the positrons (Figure 8). Considerable 
Information on the annihilation site of the positrons should become available 
from the comparison of observable time dependences of the 0.511 MeV Tine with 
theory. No such comparison has yet been done with the SMM data. 

Because of the delayed nature of both the 2.223 and 0.511 MeV lines, 
information on the timing of the acceleration of protons and nuclei can be 
best obtained from prompt nuclear deexcitation lines. The comparison of the 
time histories of such lines with those of hard X-rays of various energies can 
give information on the possible existence of multiple acceleration stages of 
energetic particles In solar flares. Because the 4 to 7 MeV channel is 
dominated by nuclear radiation (Section III A), the flux in this channel (the 
"main channel window" In the SMM data) is an excellent diagnostic of the 
timing of acceleration of the nucleonic component In flares. 

The time history of the "main channel" flux was observed for the June 7, 
1980 flare (Chupp 1982). In particular, several gamma-ray spikes were seen In 
good correlation with the seven hard X-ray spikes (Kiplinger et al. 1982). 

From the comparison of these two time histories it follows that, at least for 
the June 7 flare, the <100 keV electrons and >10 MeV nuclei were accelerated 
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In very clos^. Mme proximity (less than a few seconds). There Is, 
nevertheless, a delay of approximately 2 sec between the hard X-ray and 
gamma-ray peaks. Bal (1982) suggests that this delay Is due to two-step 
acceleration, where the first step accelerates the <100 keV electrons and the 
second step accelerates the mildly relativistic electrons and the nuclei (see 
also Bai and Ramaty 1979). On the other hand, Chupp (1982) attributes the 
delay to the difference in propagation time of >10 MeV/nucleon nuclei and <100 
keV electrons along a Mgnetic arch of length ~ lO^^cm. In any case, the 
rapid decay I'- 2 seconds) of the main channel emission for the June 7, 1980 
flare requires a sufficiently high ambient density (>10l3cm“3) and this 
provides additional support for the validHy of the thick- target model for 
this flare. 

The time dependences of the various hard photon emissions for the June 7, 
1980 flare were quite different from those of the longer duration August 4, 
1972 flare. For the latter, the X-ray continuum above 350 keV reached peak 
strength a few minutes later than the continuum above 30 keV and the 2.22 MeV 
line profile was better explained when the neutron-production time profile was 
assumed to be similar to the time profile of the >350 keV continuum rather 
than the <100 keV continuum (Bai and Ramaty 1976c Bai 1982). These results, 
together with earlier X-ray observations (Frort and Dennis 1971) and gamma-ray 
measurements (Hudson et a1. 1980, Willet et al. 1982) demonstrate that not all 
energetic particle populations in flares are accelerated at the same time. 

The present status of the existence of multiple acceleration phases in flares 
has been reviewed by Bai (1982). 


C> The Photospheric ^He Abun^nce 

T^e abundance of ^He tn the solar atmosphere Is of considerable 
astrophy steal interest. Along with and ^He, ^He Is formed by 
nucleosynthesis In the big bang (Wagoner 1973). In addition, stellar 
evolution should Increase the ^He abundance# In particular, the Sun should 
have burned Into ^He any amount of deuterium ft originally had but the ^He 
should not have been further burned Into (Geiss and Reeves 1972 and 
references therein). Thus, a measurement of the solar ^He abundance provides 
an upper limit on the protosolar abundance, which, In turn, provides 
Information on nucleosynthesis In the big bang and on whether the universe Is 
open or closed (Gott et al. 1974). 

%e has been observed In the solar wind where the ^He/^He ratio Is of the 
order of a few times lO"^ (Geiss and Reeves 1972). Hall (1975) determined 
spectroscopically a ^He/^He ratio of (4 ± 2)xl0"^ In a solar prominence. 

There is, however, no direct observation of ^He In the photosphere. The solar 
gamma-ray observations can set limits on the photospheric ^He abundance. This 
can be done In two ways. First, for a ^He/H ratio much larger than 5x10"^, 
the 2.22 MeV line fluence would be much lower relative to the 4 to 7 MeV 
fluence than observed. Second, the ^He/H ratio must have an upper bound at a 
value not much higher than BxlO"^ because otherwise the delay between neutron 
production and 2.223 MeV photon release would be shorter than observed. There 
are as yet no firm values on the upper limits on the ^He/H ratio from the SMM 
data, but a safe preliminary limit would be ^He/H < 2x10”^. 

In addition to setting an absolute upper bound on the photospheric ^He/H 
ratio, the gamma-ray data also limit any possible variability of this ratio In 
time and with position on the solar surface. From the small variability of 
<|)(4-7)/(j>(2.22) from flare to flare (Table 2), it follows that ^He/H should be 
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constant to better than a factor of 2. 

As discussed in the Introduction, very high ^He/^He ratio have been 
observed in solar energetic particles {see Ramaty et a1. 1980 for a review of 
these observations). It is now believed that these enhancements are not of 
nuclear origin, but result from selective heating and acceleration (Fisk 
1978). Finite ^H/^H ratios in energetic solar particles averaged over several 
solar flares have been presented (Anglin, Dietrich and Simpson 1973b, Hurford, 
Stone and Vogt 1975). But the very large uncertainties in these measurements 
and the possibility of instrumental contamination preclude a definite 
conclusion regarding the positive detection of secondary nuclear products in 
energetic solar particles. It appears, nonetheless, that the bulk of the 
nuclear reactions are produced by the flare accelerated particles that remain 
trapped at the Sun. The particles that escape from the Sun and are observed 
in the interplanetary medium are devoid of any measurable amount of nuclear 
spallation products. 

D. Beaming of the Energetic Particles 

As we have seen in Section II D, gamma- ray line observations can give 
unique information on the beaming of the energetic particles. Shifts in the 
energies of narrow lines are indicative of such beaming (Ramaty and Crannell 
1976), but these effects are probably difficult to measure with low resolution 
spectrometers. Another effect of beaming is the narrowing of the broad 
lines. These lines, produced by energetic heavy nuclei, are Doppler broadened 
by both the /elocity spread and the angular distribution of the particles. In 
the case of a beam, however, the latter effect is greatly reduced and hence 
broad lines (Figure 9) can mimic narrow lines. For the same reason, the ^Be 
and \i lines shown in Figure 13 are much narrower if produced by a-particles 
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1n a beam than by such particles with an Isotropic distribution. 

The energies of these lines would then provide direct Information on the 
angle between the beam and the direction of observation. A flare model tn 
which the gamma**ray lines would be produced by a beam of energetic particles 
Is that of Colgate (1978). Future gamma-ray line observations and more 
refined analysis of the SMM data should produce much new Information on 
energetic particles beams In solar flares and hence on the flare model. 


IV, SUMMARY 

Gamma-ray lines are the most direct probe of nuclear processes In the 
solar atmosphere. The line observations from a number of flares, made with 
spectrometers on OSQ-7, HEAO-1, HEAQ-3 and SMM, are consistent with reactions 
produced by flare accelerated particles of energies greater than several 
MeV/nucleon. These reactions involve the production of neutrons, emitters, 
Tf mesons and excited nuclear levels, all of which lead to observable gamma-ray 
line emission. 

The solar gamma-ray line observations can give information on the timing 
of the nucleonic component in flares, through measurements of the light curves 
of prompt lines, on the energy spectrum, number and energy content of these 
particles, through line ratios and line fluences, on the site of the nuclear 
reactions, through selective attenuation of lines from limb flares and the 
spectrum and time dependence of e'^-e’ annihilation radiation, on the geometry 
of particle beams, through line shapes and Doppler shifts, on the photospheric 
^He abundance, through the time dependence and fluence of the 2.223 MeV line, 
and on chemical compositions of both the ambient medium and the energetic 
particles, through gamma-ray line ratios. 

At the time of this writing (early 1982), only limited portions of the 
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SMM gamma-ray data are available for analysis and hence only few hard 
conclusions can be drawn from them (see also Chupp 1982). 

The acceleration of protons and nuclei to energies above an MeV, In at 
least some flares, takes place In a time Interval less than a few seconds. 

This sets Important constraints on flare acceleration mechanisms which have 
not yet been fully explored. For at least some flares, previous Ideas on two 
phases of acceleration, which Involve long delays (>1 minute) between the 
acceleration of MeV nuclei and X-ray producing electrons, are not valid. But 
there Is evidence that the > 10 MeV protons are accelerated later than the < 
100 keV electrons. 

The 2.223 MeV line Is strongly attenuated in limb flares. This provides 
direct observational confirmation for neutron capture In the photosphere. 
Further confirmation of this process comes from the precise measurement of the 
energy of this line (with the high resolution Ge detector on HEAO-3) and from 
the observed delay between the 2.223 MeV flux and the flux of prompt nuclear 
radiation. The 2,223 MeV line observations also Indicate that the 
photospheric ^He abundance is about 5xl0“5 relative to M by number, and that 
It does not vary much with position on the Sun. 

The gamma-ray line emission Is produced in thick-target interactions, 

I.e. by energetic protons and nuclei which slow down in the solar 
atmosphere. The absence of nuclear fragments (^h, ^h, l 1. Be, B) in the 
fluxes of the interplanetary particles Indicates that thin- target interactions 
do not produce many gamma rays. This Is probably the reason for the lack of 
correlation between the number of particles responsible for gamma-ray line 
production and the number observed in the interplanetary medium. 

The ratio of the 4 to 7 MeV photons fluence to the fluence In the 2.223 
MeV line Is a strong function of the energy spectrum of the accelerated 
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particles. This spectrum, as deduced from the gamma-ray observations » does 
not vary much from one flare to another, and is, within rather broad error 
ranges, similar to the particle energy spectra observed in the interplanetray 
medium. This argues for the same acceleration mechanism for both the 
gamma-ray producing particles and the interplanetary particles. In this case, 
the lack of correlation between the absolute numbers of particles in the two 
populations could be due to the varying escape conditions of energetic 
particles from solar flares. Multiple acceleration mechanisms, however, 
cannot be ruled out at the present time. 

Much additional information on solar flares and on energetic particle 
acceleration therein is expected from the detailed analysis of the already 
available SMM data and from new data anticipated from the SMM spectrometer as 
well as from other spectrometers that hopefully will be flown during the next 
solar maximum toward the end of the 1980's. Of particular interest would be 
the observation of solar gamma-ray lines with high spectral resolution which 
could provide unique information on such questions as the beaming of the 
energetic particles, the temperature of the energetic particle interaction 
site and the compositions of the ambient medium and the energetic particles. 
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TABLE 1 


Elemental and Isotopic Abundances 


Isotope 

Ambient 

Energetic 

Particles 

Particles 


1. 

1. 

fUe 

0.07 

4.15 X 10“J 

0.15 


1.07 X lO'^^ 


4.64 X 10"® 

1.28 X 10"® 


9.0 X 10"® 
3.46 X 10"J 

2.14 X 10"i 


8.57 X 10"' 


6.92 X 10“J 

2.14 X 10“3 

IBq 

1.38 X 10“® 

4.28 X 10"® 

f^Ne 

9.0 X 10"® 

2.14 X 10"J 


1.0 X 10"®^ 

2.57 X 10"® 


2.28 X 10"® 
3.11 X 10"® 

4.28 X 10"® 

IS 

6.42 X 10"J 

ofwg 

4.01 X 10"® 

8.14 X iO"® 

?SMg 

4.43 X 10"® 

8.49 X 10"® 

27a1 

3.18 X 10"® 

5.35 X 10"® 

??Si 

3.46 X 10-5 

6.42 X 10"^ 

?%i 

1.80 X 10"® 

3.21 X 10"® 

?9si 

1.18 X 10"® 

2.14 X 10"® 


1.80 X 10"® 

1.07 X 10-J 


7.61 X 10"' 

4.71 X 10"® 

onAr 

3.39 X 10"® 
6.23 X 10"' 

2.14 X 10"® 


4.28 X 10"® 

c?Ca 

2.28 X 10"® 

4.28 X 10"® 

®JCr 

4.15 X 10"7 

2.14 X 10"® 

IS 

1.94 X 10"® 

6.85 X 10"® 

56pe 

3.11 X 10"® 
7.61 X 10"' 

1.07 X irj"^ 

IS 

2.57 X 10"® 

f®N1 

1.25 X 10"® 

2.14 X 10"® 

60^1 

4.84 X 10"' 

8.57 X 10"® 


Note: Two sets of energetic particle abundances are discussed in the text; 
tPl - Energetic particle abundances given in Column 3. 

EP2 - Energetic particles with the same abundances as the ambient 
medium, Column 2. 
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FI are 

TABLE 2 

Gamma-Ray Flares 
Fluences, l)(pfjotons/cm^) 

2.22 MeV 4-7 MeV 4.44 MeV 6.13 MeV 0.511 MeV 

4>(4-7) 

J(2‘.2Z) 

Location 

1972, 

Aug 4 

155±12^ 

105±11^ 

17±5^ 

17±5^ 35±ll' 

0.68±0.09 

EOS 

N14 

1972, 

Aug 7 

No data available during 

time of maxi mum emission 

— 

W38 

N15 

1978, 

July 11 

240±70^ 

171^ 

43±17^ 

— 

0.71 

E43 

N18 

1979, 

Nov 9 

38±9® 

50±4® 

— 

— 

1.32i0.33 

EOO 

$16 

1980, 

June 7 

6.6±1® 

11.5±0.5® 

— - 

<2 

1.74i0.27 

W74 

N12 

1980, 

June 21 

No data available during 

time of maximum emission* 

— 

W91 

N17 

1980, 

July 1 

3.3i0.5® 

3.1±0.4^ 

— 

--- 0.9±0.4 

0.94±0.19 

W37 

S12 

1980, 

Nov 6 

10.3±1.3^ 

14.8±0.8® 

— 

<2 

1.44±0.2 

E74 

S12 

1981, 

Apr 10 

13.5±1® 

18.6±1.6® 


<6.6 

1.38±0.16 

W37 

N09 

1981, 

Apr 27 

11.7±2^ 

118±2® 

— 

— — 

10.U1.7 

W90 

N16 


1. Chupp {1975) 

2. Ramaty, Kozlovsky and Suri (1977) 

3. Hudson et aV. (1980) 

4. Theoretical 

5. T. Prince (private communication 1981) 

6. Chupp (1981) 

*Delayed 2.22 MeV and 0.511 MeV line, and high-energy neutrons were observed from this flare 
(Share et al. 1980, Chupp and Forrest 1981). 



TABLE 3 


Spectral Parameters. Total Energies and Number of Particles of Ganna Ray Flares 


FI are 


Q(4-7) 

^n 

aT 

W(>lMeV) 

(erg) 

Np(>10 Mev) 

N«c!p 

1972 Aug 4 

1 

0.16 

0.019 

2.5)^1030 

1.3xl03‘^ 

3x1035 

1978 July 11 

0,9 

0.15 

0.020 

1.8x1030 

1.0xl03‘+ 

— 

1979 Nov 9 

1 

0.30 

0.014 

9x1023 

3.4x1033 

... 

1980 June 7 

0.6 

0.24 

0.015 

2x1023 

8.5x1032 

1031 

1980 July 1 

0.95 

0.21 

0.016 

5x1028 

2.3x1032 

«... 

1980 Nov 6 

0,6 

0.20 

0.017 

2x1029 

1.0x1033 

— 

1981 Apr 10 

0.95 

0.30 

0.014 

3.5x1029 

1.3x1033 

— 

1981 Apr 27 

~ 0 

- — — 

0.019 

1.4x1030 

7.3x1033 





0.511 


Flare 

Duration 

aT 

Q+/Q(4-7) 

Aug 4. 1972 

553 sec 

0.019 

0.46 

June 7, 1980 

50 sec 

0.015 

0.41 

July 1, 1980 

60 sec 

0.016 

0.41 


Apr 27, 1981 2000 sec 


0.019 


0.46 
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TABLE 4 

MeV Line Fluences 

1>(0.5l)(Photons/cin2) 
Calculated Observed 


48 

^0.51 

35±11 

4.7 

^0.51 

<2 

1.3 

^0.51 

0.940.4 

54 

^0.51 

<10 
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Figure Captions 

1. Interplanetary protons and a-part1cles observed on June 7, 1980 (R. 
HcGuIre, private communication 1980). The curves ar-e from equation (11) 
with the same aT for protons and a-par*t1cles. 

2. Neutron yields In the thln-tar-get and thick-target models calculated with 
the energetic particle spectra of equation (10) (exponential In rigidity) 
and equation (11) (Bessel function). The ambient medium and energetic 
particle compositions ar*e from Table 1. Ng 5 g^p(>E) and Rp(>E) are, 
respectively, the numbers of protons of energies greater than E that escape 
from the thin-target region or are Incident on the th1ck-tar*get, n^^T Is 
the product of the hydrogen density In the thin-target region and the 
particle escape time from this r-eglon. In the thick-target calculations 
the particle energy loss rate Is given by equation (7) appropriate for for 
a neutral medium. 

3. Partial neutron production r*ates in the thick-target model for Bessel 
function energetic particle spectra. The various curves give the neutron 
production by energetic protons (P), a-part1cles (a), C, N and 0 nuclei 
(CNO) and nuclei from Ne through Fe (Ne-Fe) Interacting with the ambient 
medium. All the other parameters are as In Figure 2. 

4. Energy deposition per neutron produced 1n a thick-target model for Bessel 
function energetic particle spectra. The various curves give the energy 
deposited by particles of energies greater than the indicated values. The 
compositions of the ambient medium and ener-getlc particles are given in 
Table 1 and the energetic particles are slowing down in a neutral medium. 

5. Probabilities for neutron escape, decay, capture on protons and loss on 
^He in the solar atmosphere (solid lines); and photon yields per neutron 
(dashed lines). The parameter 9 Is the angle between the Earth-Sun line 




and the heliocentric radial direction through the flare. The ratio ^He/ll 
is the photospheric ^He abundance, and is the energy of the neutrons. 

The initial neutrons are assumed to be released isotropically above the 
photosphere (from Wang and Ramaty 1974). 

6. The ratio of the positron yield, Q+, to the neutron yield, Q^j, for Bessel 
function energetic particle spectra (equation 11). The composition of the 
ambient medium and energetic particles are given in Table 1 and the 
energetic particles are slowing down in a neutral medium. 

7, The instantaneous fractional positron production rate and the time 
integrated fractional positron yield in a thin- target model for a burst of 
e"** emitter production (6 function) at t = 0, The energetic particles 
spectrum is given by equation (9) and the compositions are close to those 
of Table 1. The results of this figure do not depend strongly on the 
interaction model used. 

8, The slowing down time, tg, plus positronium formation time, tpos» of 

positrons of initial energies E+ in an ambient medium of temperature lO^K, 
degree of ionization n and density 10^1 or 10^3 ^^,-3, unear parts 

of the curves, tg » tp^gj tg is essentially independent of temperature and 
does not depend much on n; fpos» however, depends strongly on these 
parameters (see Bussard, Ramaty and Drachman 1979 for more details). The 
calculations of Figure 8 have been carried out by J. M. McKinley (private 
communication 1981). 

9. Prompt nuclear gamma-ray spectrum from the interactions of energetic 
particles with the solar atmosphere in a thin-tar-get model. The 
composition of the ambient medium is given in Table 1. The energy spectrum 
of the energetic particles is given by equation (9) and their composition 
is the same as that of the ambient medium (EP2). Not shown in this figure 
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are the delayed Tines, at 2.223 MeV from neutron capture and at 0.5U MeV 
from positron annihilation. In the August 4, 1972 flare, these lines were 
10 and 2 times more Intense, r’espectively, than the vine at4.44 
MeV. (From Ramaty et al. 1980). 

10. The ratio of the narrow 4.44 MeV line yield to the neutron yield for 
energetic particle spectra given by equation (11). The compositions of the 
ambient medium and ener*get1c particles are given In Table 1. 

11. The ratio of the photon yield In the 4 to 7 MeV channel to the neutr'on 
yield for energetic particle spectra given by equation (11). The 
composition of the ambient medium Is given In Table 1. The curves EPl ar-e 
for energetic particles of composition as given In Table 1 while the curves 
EP2 are for an energetic particle composition which Is the same as that of 
the ambient medium In Table 1. 

12. The r*at1o of the narrow 4.44 MeV line yield to the 4 to 7 MeV channel 
yield. The rest of the par’ameters are the same as In Figure 11. 

13. Photon spectrum of the prompt \l and ^Be lines for an a-part1c1e beam 
confined to a cone with half opening angle 5* or 20^*, and for an Isotopic 
distribution; is the angle between the beam and the direction of 
observation. (From Kozlovsky and Ramaty 1977). 
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